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abstract 

In this paper we shall study quantum critical behavior of lattice model of unconventional supercon- 
ductors (SC) that was proposed in the previous papers. In this model, the Cooper-pair (CP) field is 
defined on lattice links in order to describe c?-wave SC. The CP field can be regarded as a U(l) lattice 
gauge field, and the SC phase transition takes place as a result of the phase coherence of the CP field. 
Effects of the long-range Coulomb interactions between the CP's and fluctuations of the electromag- 
netic field are taken into account. We investigate the phase structure of the model and the critical 
behavior by means of the Monte Carlo simulations. We find that the parameter, which controls the 
fluxes (vortices) of the CP, strongly influences the phase structure. In three-dimensional case, the 
model has rich phase structure. In particular there is a "monopole proliferation" phase transition 
besides the SC phase transition. Depending on the parameters, this transition exists within the SC 
phase or takes place simultaneously with the SC transition. This new type of transition is relevant 
for unconventional SC's with strong spatial three-dimensionality and to be observed by experiments. 
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1 Introduction 



Sonic of the experiments on the high-T^ materials indicate that the superconducting (SC) pliase 
transition at vanishing temperature (T) is a quantum phase transition [1] and it belongs to the 
universality class of the four-dimensional (4D) XY model[2]. On the other hand for constructing a 
Ginzburg-Landau (GL) theory of the (i-wave SC phase transition, the Cooper-pair (CP) field must 
be put on lattice links instead of site because its condensation has different signature depending on 
the direction of the CP (see later discussion). It is also suggested that the SC transition at T = is 
a phase decoherent phase transition and the SC transition of the liigh-T'c cuprates is more like the 
Bosc-Einstcin condensation than the BCS transition for the coherent length of the high-Tg SC's is 
very short compared to the ordinary metallic SC's. 

In the previous papers[3, 4, 5] from the above point of view, we have introduced a CP field Vx^i 

Vx,i ~ {■tlJa:,1'4)x+i,l - V'xaV'x-Fi.r) (1-1) 

where ipx,a{'^ =Tj i) is the electron at site x with spin a and i = 1, • • • , d is the spatial direction 
index (it also denotes the unite vector). As we study a quantum phase transition at T = 0, the 
GL theory is defined on the d + 1 (d-spatial and 1-time) dimensional hypercubic lattice. The CP 
field Vx,i has electric charges at x and x + i and it also interacts with the transverse electromagnetic 
(EM) vector potential Ax^i- 

For the ordinary s-wave SC, the GL theory defined on lattice plays an important role for study 
on the critical phenomena. There the Cooper pair (px ~ {i^x-fi^xi " i^xi^'x^) is put on lattice site, 
and its condensation induces the SC phase. This GL theory is called gauge-Higgs model in the 
elementary particle physics and plays an important role in the unified theories. Recently the gauge- 
Higgs models have renewed interests and have been studied intensively. In particular, dynamics 
of multi-flavor Higgs models and their phase structure, order of Higgs phase transitions, eff'ect of 
the Berry's phase term, etc. are important topics. (See Refs.[6, 7].) These studies are related to 
multi-gap SC's, deconfined critical points in strongly-correlated electron systems, superfluid and SC 
in liquid hydrogen at low-T and in high pressure, etc. We expect that the gauge-Higgs model with 
the link Higgs field Vx^i also plays an important role in various fields of physics, and to clarify its 
dynamics gives useful insight to critical phenomena. In the present paper, we focus on its application 
to the unconventional (UC) SC. 

In the previous paper [5], the following action of the GL theory of the UCSC was derived from a 
microscopic Hamiltonian by using the path-integral methods (see Fig.l), 

5gl = ^[ci E Plu{x)-C2 Y.{Ux,^V:^^^iUx+i,^Vx,i + C.C.) 
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-d2 ^ (C/...C/.+.,K*+,,14j + C.C.) - C3 ^ {V.,^V:+,^,V,+,,^V:^, + C.C.) 



(1.2) 



where F^y{x) = A^j^^ + A^+^^y - A^^u - A^+i,^^ v = 0,1, ...,d) and Ux^p, = e*-^"'*". The foUowing 
comments on ^gl are in order. 

1. The 0-th component vector potential A^fl mediates the long-range Coulomb interaction be- 
tween the CP's. 

2. S'gl is invariant under the following gauge transformation, 



(1.3) 



where ax is an arbitrary function of x. 

3. ci ~ /3 in ^GL are parameters and the overall factor g plays a role of 1/H and controls quantum 
fluctuations[2, 5]. 

4. The C2-term represents the hopping of Vxa, whereas the d2-term determines the relative phase 
of the adjacent Vx/s, e.g. the negative d2 enhances d-wave condensation of Vx,i- 

5. The C3-term in ^gl controls plaquette-flux (i.e., plaquette- vortex) of Vx^- 

6. For large C3, we expect that monopole configurations of Vx^ are suppressed and effects of the 
compactness of Vx,i becomes negligibly small. 

7. Parameters a and (3 control the mean value of Vx,i and its fluctuations. 
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Figure 1: 5gl in Eq.(1.2). 
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In the previous paper[5], we investigated the phase structure of the models (1.2) for ci = C2 = 
—d2 = C3 = 1 as varying the parameter g and found existence of second-order SC quantum phase 
transitions. In the = 3 cases, however, we found another phase transition within the SC phase 
besides the SC transition and we discussed its physical meanings. Though the high- Tc cuprates 
are quasi-two-dimensional and have layered structure, critical behavior is governed by its three 
dimensionality. Therefore the above additional phase transition might be observed in experiments. 

In this paper we shall continue the previous study and investigate phase structure and quantum 
critical behavior of the model for various parameters regions. In particular, we shall discuss and 
verify that the locations of the two second-order phase transitions get closer with each other and 
finally they merge to a single first-order phase transition as the value of C3 is increased. 

This paper is organized as follows. In Sec. 2, we shall study the London limit (the CP field 
\Vx^i\ = 1) of the model in (3 + 1)D. For the case ci = C2 = — ^2 = C3 = 1, there are two second-order 
phase transitions as observed in the case investigated in the previous paper. We discuss the physical 
meanings of the phase transitions. Then we study how the phase structure will change as the value 
C3 is increased. In Sec. 3, we shall study the case a = (3 = 5 and show that various phase structures 
appear as the value of C3 is varied. There new type of phase transitions appear, which are relevant 
for SC materials with strong three-dimensionality and to be observed by experiments. Section 4 is 
devoted for conclusion. 

2 London limit in (3+l)D 

In the previous paper [5], we reported the results of study on the a = p = 5 and a = /3 = 10 
cases. We investigated the case with the parameters ci = C2 = — ^2 =03 = 1, and found that in 
(3-|-l)D there are two phase transitions at Qc and g[. {qc < g'c)- By the calculation of the gauge-boson 
mass Mq, we concluded that the first transition at Qc is the genuine SC phase transition, which is 
connected with the proliferation of vortices of the CP field Vx^i, whereas the second one at g'^ is 
related to the monopoles proliferation of T4,i. Furthermore we inferred that the locations of the two 
phase transitions would be getting closer and finally they would merge into a single first- order phase 
transition for larger values of C3 > 1. Recently this kind of phenomena of the phase transition were 
observed in simpler gauge-theory systems [7]. 

In this section we shall study the system S'gl in the London limit \Vx,i\ = 1. Then the system 
given by Eq.(1.2) can be regarded as a U(1)nc xU(l)c gauge model where the EM vector potential 
corresponds to the (3 + 1)D noncompact U(l) gauge field, whereas the CP field Vx,i is the 
3D compact U(l) gauge field. We shall study the system and determine its phase structure by 
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calculating "internal energy" E, "specific heat" C, the magnetic penetration depth, etc. E and C 
are defined as follows, 

E=-{Sgl)/L^ C = ((5gl - (SGL))')/L^ (2.1) 

where is the system size. The gauge-boson mass Mq is nothing but the inverse magnetic pene- 
tration depth and is obtained by calculating the correlation function of the magnetic field Fij {x) [8] . 
Other interesting quantities are monopole density of Ux,i and Vx,i, pu and pv, which measure the 
strength of the topologically nontrivial fluctuations of the link gauge fields[9, 3]. pu must have 
strong correlation to Mq and have small values, but sometimes pu gives clearer signal of the phase 
transition (see later discussion). As the CP V^^i is defined on lattice link in the present system, vor- 
tices (more precisely vortex lines) in the SC phase can be terminate at (anti)monopoles. Therefore 
pv is a very specific quantity in the present model Sgl and does not exists in other systems like 
the XY model. However as C3 is getting large, Vx^i tends to have ^^pure-gauge configurations" like 
Vx,i ~ e'^"'+*e*^== [10] and pv decreases very rapidly. 

We first study the case ci = C2 = — ^2 = C3 = 1. The calculations of C and the inverse penetration 
depth Mq are shown in Fig. 2. From the result of C, it is obvious that there are two second-order 
phase transitions at Qc ^ 0.37 and g'^ ~ 0.405, since the each peak develops as the system size L is 
getting large. From the calculation of Mq, we can see that the first transition at gc is the genuine 
SC transition. In order to see the physical meaning of the second transition at g'^, we measured the 
U and V monopole densities pu, pv- See Fig.3. pu has small values for Ux,ix is the noncompact 
gauge field and it changes the behavior at both g^ and g'^. pv has finite (fairly large) values in the 
region between gc and g'^ and it starts to decrease at g'^. This indicates that the transition at g'^ 
corresponds to the monopole proliferation-suppression of the CP field Vx^i, i.e., though the density 
of vortices decreases at the first phase transition point gc, short vortex lines, which terminate at 
(anti)monopoles, survive in the region between gc and 5^. 

Let us turn to the case C3 > 1. As explained in the above, effects of the compactness of the 
"gauge field" Vx^i is suppressed for large C3, and the pure-gauge configurations Vx,i ~ e'^^^+^e'^* 
dominate the path integral. Then one may naively expect that the system is getting close to the 
4D XY model. On the other hand, the dynamics of Vx^i is getting close to that of the noncompact 
U(l) gauge field like Ux,i because of the suppression of the ^-monopoles, and as a result the system 
becomes U — V symmetric. Then one may expect that two second-order phase transitions in the 
C3 = 1 case tend to get close with each other and finally they merge into a single phase transition. 
From the above point of view, it is very interesting to investigate the cases of larger value of C3. 

We first show the results of the C3 = 2 case. Internal energy E is shown in Fig. 4. These results 
show that there exists a finite jump in E &i g = 0.355 and its discontinuity gets sharper as the 
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Figure 2: (color online) C for the London limit with C3 = 1. There are two peaks both of which 
have system-size dependence (left). Gauge-boson mass Mq. It develops nonvanishing value at the 
first critical coupling gc ~ 0.37 (right). 
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Figure 3: Monopole densities in the London limit with C3 = 1. pu (left) and pv (right). In the 
region between = 0.37 and g'^ = 0.405, pv has finite values. 
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system size L is increased. This means that there exists a first-order phase transition at Qc — 0.355. 
No other anomalous behavior is observed in E and C as g is increased further. 




Figure 4: (color online) Internal energy E for C3 = 2 in the London limit. As the system size is 
getting larger, the discontinuity at gc = 0.355 becomes sharper. 

We expect that two second-order phase transitions, which exist for C3 = 1, merge into the first- 
order phase transition in the present C3 = 2 case. In order to verify the above expectation, we 
measured the monopole densities pu and pv- See Figs. 5 and 6. It is obvious that both pjj and py 
exhibit discontinuity at the phase transition point gc = 0.355 and are vanishingly small for g > gc- 
We also measured the gauge-boson mass Mq. The result is shown in Fig. 7. Mq does not show sharp 
discontinuity at the critical point, but it starts to increase from zero at gc — 0.355. 
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Figure 5: (color online) The U-monopole density pjj for C3 = 2 in the London limit. As the system 
size is getting larger, the discontinuity at gc — 0.355 becomes sharper. 

We have also investigated the case of the London limit with C3 = 4. We show the results in Fig. 8. 
There is a first-order phase transition at 5 = 0.320 and behavior of i?, pu and pv have similar 
behavior with those of the C3 = 2. However Mq shows a sharp discontinuity at the critical point in 
the present case. 
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Figure 6: (color online) The V-monopole density pv for C3 = 2 in the London limit. As the system 
size is getting larger, the discontinuity at gc — 0.355 becomes sharper. 
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Figure 7: The gauge-boson mass Mq for C3 = 2 in the London limit. 
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Figure 8: (color online) Internal energy E and gauge-boson mass Mq for C3 = 4 in the London limit. 
Both quantities exhibit sharp discontinuity at critical point gc — 0.320. 
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3 Phase structure of a = (3 = 5 case 



In the previous section, we showed that the two second-order phase transitions merge into a single 
first-order transition as C3 is increased. This behavior is expected from the recent studies on more 
tractable gauge-Higgs models. In this section, we shall study the system with a — f3 = 5 and see 
how the phase structure changes as the parameter C3 is varied. In this case Vx^i is parameterized 
as Vx,i = r^^iC^^'''^ and the amplitude rx,i is also a dynamical variable and fluctuates around its 
mean value. This may change the critical behavior of the model. In fact, it is known that the phase 
structure of the 3D U(l) Higgs model depends on the potential term of r^j.^H]- 

The system with C2 — ~ c^i — I has a similar phase structure to that of the London limit 
with C2 = —d2 = C3 = 1, i.e., it has two second-order phase transition points at gc ~ 0.665 and 
g'^ — 0.778. The detailed results have been given in Ref.[5]. See Figs. 9 and 10 for C, r — (vx.i) and 
the gauge-boson mass Mq. The superfluid density r = {rx.i) changes its behavior at both critical 
couplings. The mean-field approximation predicts (r — rg) cx Mq, where tq — (rx.i) for g < gc- 
The data in Fig. 10 shows the small but finite deviation from (r — tq) cx Mq due to the quantum 
fluctuations of Vxa and Ax^i- 




Figure 9: (color online) Two second-order phase transitions for C3 = 1, a — f3 ^ 5. 



In the rest of this section, we shall consider the cases with C3 = 2 and C3 = 4. We first show the 
internal energy and AIq for the case of C2 = —d2 ~ 1 and C3 = 2 in Fig. 11. Calculation of E indicates 
that there exist two first-order phase transitions at gc = 0.63 and g'^ — 0.65. We have not observed 
this kind of behavior for the system in the London limit. Mq starts to develop at gc and therefore 
the first transition at gc is the SC phase transition. In order to verify the above conclusion that two 
first-order phase transitions exist, we measured pu and pv- See Fig. 12. The both quantities exhibit 
sharp discontinuity at the critical couplings gc and obtained by E. However, the discontinuity 
of Pu (pv) at gc {g'c) is larger than that at g'^ (gc)- In the region between gc and g'c, pv has finite 
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Figure 10: (color online) r — {r^.i) for C3 = 1, a = /3 = 5(left). Mq and r as a function of g(riglit). 



values as in the previous cases. 
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Figure 11: (color online) Internal energy E and gauge-boson mass Mq for a — [3 — ^ and C3 = 2. 
The data of E indicate the existence of two first-order phase transitions at gc — 0.63 and g'^ — 0.65. 
Mq develops at the first phase transition at gc = 0.63, the SC transition, whereas it shows no clear 
change in behavior at g'^ — 0.65. 



Let us turn to the C3 = 4 case. E in Fig. 13 shows that there is single first-order phase transition 
as in the London limit, r = {r^.i) and pv also show a discontinuity at gc — 0.55. In the present 
case, two first-order phase transitions at C3 = 2 merge into a single first-order phase transition. 

4 Conclusion 

In this paper we have studied the phase structure and critical behavior of the GL theory that we 
proposed earlier. In particular, we focused on the isotropic 3D case and investigated its quantum 
phase structure. In the previous paper, we found that the 2D GL system has a second-order SC 
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Figure 12: (color online) pu and pv for a = f3 = 5 and C3 = 2. The data support the existence of 
two first-order phase transitions at gc — 0.63 and g'^ — 0.65. 
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Figure 13: (color online) E, r = {r^.i) and pv for a = /3 = 5 and C3 = 4. The data support single 
first-order phase transition at gc — 0.55. 
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phase transition that is in accordance with the experiments. 

Though most of the materials of unconventional rf-wave SC's have layered structure, properties 
of the phase transitions are governed by their three-dimensionality, e.g., the finite-T phase transi- 
tion is possible by the three-dimensionality of the materials. Then we studied the 3D system and 
found that interesting phase structure appears depending on the magnitude of the parameters. For 
materials with strong three-dimensionality, the results obtained in the present paper predict one of 
the following alternative possibilities; 

1. There exists CP monopole proliferation (suppression) phase transition within the SC phase 

2. Single first-order SC phase transition exists accompanying monopole transition of the CP 

As the real material of high-Tc SC have strong anisotropy, signal of the phase transition within the 
SC phase and first-order phase transition might be too weak to be observed. However the prediction 
itself is very interesting. 

In the GL theory obtained from the canonical microscopic model of the high-Tc SC like the t-J 
model, the parameter C3 is an increasing function of the antiferromagnetic (AF) exchange coupling 
J, the amplitude of the resonating- valence-bond configuration and the hole concentration [12]. Then 
it is interesting and also important to take into account the effects of the AF background, i.e., 
the coefficients in 5gl are not simple parameters but are determined by the dynamics of the AF 
background. This problem is under study and the result will be reported in a future publication. 

Another interesting problem is the relation between the magnetic penetration depth at T = 
(the inverse gauge-boson mass at T = 0) and the SC phase transition temperature Tc. In principle, 
Tc can be calculated as a function of g in the present GL theory by means of MC simulations. 
Then the exponent a, oc (Mq)", should be compared with the experimental data of the high-T^ 
cuprates. (Experiments give a ~ ^[1].) This problem is also under study and the result will be 
published in a future publication. 
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